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Abstract 

The probability density of the resistance of a two dimensional rectangular network 
between two conducting plates is calculated. The nodes form an M by N lattice, and 
each edge has a random resistance. The Monte Carlo method is used. 

1 Introduction 

A random electrical network is a network of resistors, with the resistance of each resistor 
being a random variable. Then the overall resistance of the network between any two nodes 
is also a random variable. Its probability distribution is determined by the probability 
distributions of the individual resistors, and by the structure of the network. 

We shall calculate the probability distribution of the resistance of a two dimensional 
rectangular network between two conducting plates. The nodes of the network form an M 



by N lattice, and each edge is a resistor. The resistances of the resistors are independent 
identically distributed (i.i.d.) random variables, taking the value r < 1 with probability p, 
and the value 1 with probability 1 — p. See Figure 1. 

In the special case M — 1, the network is just a line of N resistors in series. Its resistance 

TV 

is R = r ji where r 3 - is the j-th resistance. Thus the probability distribution of R is 
just that of the sum of iV i.i.d. random variables, which is a binomial distribution with 
parameters N and p. In the present case the possible values of R and their probabilities are 



R k = kr + (N-k) } P(R k )= [ p , k = l,---,N. (1.1) 



2 Formulation for M > 1 

For M > 1, we must determine R by first solving Kirchoff's equations for the currents and 
potentials. Let i,j denote the (j + l)-st node on the horizontal line i, with j = on the 
plate at the left and j = N at the plate on the right. Let <pij be the potential at node i,j. 
Finally let r^+ijj be the resistance of the resistor from i,j to i+ 1, j, etc., and let 
be the current in that resistor. Then Kirchoff's equations can be written as 

4>i+l,j ~ 4>i,j — Iij,(i+l)j r ij,{ i + 1 )j (2-1) 

We must solve (12.11) for the currents and potentials, subject to the boundary conditions 

0*o = 0, <j) lN = I. (2.2) 
Then the total current from the plate at the left to that at the right is: 

M 

/ = X)w ( 2 - 3 ) 

1=1 

The resistance between the plates is R = J -1 , which is the random variable, the distribution 
of which we seek. We normalize it by setting R = MR/N. Then r < R < 1. 



To find the probability distribution P(R) we use the Monte Carlo method. This consists 
in generating a sample network by selecting each resistance to be r with probability p, or 1 
with probability 1—p, solving (12. ip and (12.21) . and calculating / from (I2.3p . Then we obtain 
the sample value: R = MI~ 1 /N. By repeating this process many times, we generate a large 
set of sample values. We sort them into bins, which are subdivisions of the interval (r, 1) of 
the R axis. The fraction of samples in each bin, divided by the bin length, is the empirical 
probability of R at the midpoint of the bin. The mean of the sample values is the empirical 
expected value E(R). 

The network has four parameters r, p, N, and M. In our calculations we set M = N so 
R = R. Only the three parameters r, p and N remain. 

3 Results 

We calculated P(R) and E(R) for N = 2, 10, 20, 30; r = .25, .50, .75 and p = .25, .50, .75. In 
each case we used 2000 samples, and the calculated P(R) was close to a normal distribution, 
being closer the larger the value of N. Furthermore the expected value was close to the 
power law, E(R) ~ r p , again becoming closer as N increased. Some of the calculated results 
are shown in Figures 2-5. 

In Figure 2, the calculated P(R) is shown together with a closely fitting normal distribu- 
tion for r = 1/2, p = 1/2 nd N = 10, 20, 30. In Figure 3a, the calculated E(R) is shown as 
a function of p for N = 20 and three values of r. In Figure 3b, E(R) is shown as a function 
of r for N = 20 and three values of p. In each case the function r p is also shown. 

In Figure 4a, the calculated E(R) is shown as a function of p for N = 2 and r = .5 
together with r p , and the exact result for E(R). The latter was obtained by solving (12.11) 
and (12.21) for R using Mathematica, and then obtaining E(R) from the solution. The result 
is given in the appendix. Figure 4b is similar with p = .5 and r varying. 
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Appendix 



4(1 -pfp (3 + 5r) (1 -pfp 3 (3r + 4r 2 + r 3 ) 
5 + 3r + 1 + 4r + 3r 2 



= (l-p) 5 + (l-p) 4 p + p 5 r + 



(1 -p) p A (5r 2 + 3r 3 ) (1 - p) p 2 (l + 4r + r 2 + r (1 +r)) 
+ 3r + 5r 2 + 2 + 5 r + r 2 

(l-p)V (l + 4r + r 2 + r (1 + r)) (l-p)V (l + 2r + r 2 + 2r (1+r)) 



2 + 4r + r (1 + r) 



3 + 3r + r (1+r) 



(l-p)V (l + 4r + 3r 2 ) 2(1 -p)V (l + 2r + 2r 2 + r (2 + r)) 



3 + 2r + r (2 + r) 



2 + 6r 



^ (1 -pfp 2 (1 + 2r + r (1 +r) +r (2 + r)) (1 -pfp 2 (1 + 2r + r (1 + r) + r (2 + r)) 



2 + 4r + r (1 + r) 



2 + 3r + r (2 + r) 



^ (1 -pfp 2 (3r + r (1 + r) +r (2 + r)) (1 -pfp 2 (3r + r (1 + r) +r (2 + r)) 



\3„2 



3 + 3r + r (1 + r) 



3 + 2r + r (2 + r) 



^_ (1 - pf P 3 ( r + 2 r 2 + r (1 + r) + r 2 (2 + r)) (1 - pf p 3 (r + 2 r 2 + r (1 + r) + r 2 (2 + r)) 



1 + 2r + r (1 + r) + r (2 + r) 



1 + 4r + r (1 + 2r) 



(1 - p) p 4 (5 r 2 + 3 r 3 ) (1 - pf p 3 (r + 2 r (1 + r) + r 2 (2 + r)) 



2r + 3r 2 + r (l + 2r) 



l + 2r + r (l+r)+r (l + 2r) 



(1 - pfp 3 (3r 2 + r (1 + r) + r (1 + 2r)) (1 -pfp 3 (3r 2 + r (1 + r) +r (1 + 2r)) 



1 +4r + 3r 2 



1 + 2r + r (1 + r) +r (1 + 2r) 



(1 - pfp 3 (2r + r 2 + r 2 (1 + r) +r (1 + 2r)) (1 -pfp 3 (2r + r 2 + r 2 (1 + r) + r (1 + 2r)) 



4r + 2r (1 + r) 



4r + r 2 + r (2 + r) 



+ 



(1 - pfp 3 (r 2 + r (1 + r) + r 2 (1 + r) +r (1 + 2r)) 
l + 3r + 2r(l + r) 



| {I -pfp 3 (r 2 + r (1 +r) + r 2 (1 + r) + r (l + 2r)) | (1 - p) p A (3r 2 + r 2 (1 + r) + r 2 (1 + 2r)) 



l + 2r + r (l + r)+r (2 + r) 



r + 3r 2 + 2r (1 + r) 



(1 - P) P A (3 r 2 + r 2 (1 + r) + r 2 (1 + 2 r)) (1 - p) p 4 (r 2 + 2 r 2 (1 + r) + r 2 (1 + 2 r)) 



2r + 3r 2 + r (l + 2r) 



r + 2r (l + r)+r (l + 2r) 
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Figure Captions 



1. A two dimensional network with A = 10 rows and M — 10 columns of resistors. At the 
left side is a plate with potential <f> ii0 = and at the right side a plate with <f> i>10 = 1. 

2. Simulated probability density P(R) using 2000 samples, together with a closely fitted 
normal density for r = .5, p = .5 and a) A = 10, b) A = 20, c) A = 30. 

3. Simulated mean E(R) (crosses) using 2000 samples, together with the curve r p , for 

= 20. a) r = .25, .5, .75 with p varying, b) p — .25, .5, .75 with r varying. 

4. Simulated mean E(R) (dotted line), r v (thin line) and exact E(R) (heavy line) for 
N = 2. a) r = .5 with p varying, b) p — .5 with r varying. 

5. Simulated probability density P(R) as a function of R and A with r = .5 and p = .5. 



6 



Figure 1: 
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(a) 



Simulation vs Gaussian, N=10 
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(b) 



Simulation vs Gaussian, N=20 
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(c) 



Simulation vs Gaussian, N=30 
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Figure 2: 
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E[R] for varied p, r — .25, .5, .75 E[R] for varied r, p — .25, .5, .75 




Figure 3: 




Figure 4: 
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Figure 5: 
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